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COMPUTER  PROGRAM  THAT  PROCESSES  ENVIRONMENTAL 
DATA  TO  FACILITATE  CONTOURING 


1.0  IDENTIFICATION 

1.1  Title 

A program  to  find  the  high  and  low  values  and/or  find  or  interpolate  all  values 
which  are  multiples  of  a specified  value. 

1.2  Identification  Name 
HIGHLOW. 

1.3  Classification  Code 
None. 

1.4  RCC  Identification  Number 
None. 

1.5  Entry  Points 
HIGHLOW. 

1.6  Programming  Language 

Language:  CDC  3600/3800  Fortran. 

Routine  type:  program. 

Operating  System:  Drum  Scope  2.1. 

1.7  Computer  and  Configuration 
CDC-3800. 

1.8  Contributor  or  Programmer 

Marilyn  L.  Blodgett,  Code  8122MB,  Applied  Ocean  Acoustics  Branch,  Acoustics 
Division,  written  for  the  Environmental  Sciences  Group,  Acoustics  Division. 

1.9  Contributing  Organization 

NRL  — Naval  Research  Laboratory,  Washington,  DC  20375. 

Manuscript  submitted  December  27,  1978. 
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1.10  Program  Availability 

If  supplied  with  a magnetic  tape,  the  Environmental  Sciences  Group,  Acoustics 
Division,  will  make  a copy  of  this  program  available. 

1.11  Verification 


This  program  has  been  used  and  tested  by  the  Environmental  Sciences  Group, 
Acoustics  Division. 

1.12  Date 

June  1978. 


2.0  PURPOSE 

2.1  Description  of  the  Routine 

This  program  reads  the  bathymetric  or  magnetics,  lata  collected  by  an  oceano- 
graphic or  geophysical  experiment  from  a r .agnetic  tape  or  tapes  and  finds  the 
high  and  low  values  and/or  all  readings  which  are  multiples  of  a specified  value 
(such  as  every  100  meters  or  every  100  gammas).  It  also  checks  to  see  if  there 
are  multiples  of  the  specified  value  between  any  two  readings.  If  so,  the  program 
will  interpolate  the  required  fixes.  The  input  tape(s)  containing  the  bathymetric 
and/or  magnetics  data  is  written  in  the  format  recommended  by  the  National 
Academy  of  Sciences.  There  is  one  logical  record  (of  80  characters)  for  each  data 
point.  The  input  tape(s)  may  also  contain  navigational  data,  but  each  data  type 
is  separated  by  an  end-of-file  mark  with  a double  end-of-file  mark  at  the  end  of 
all  the  data. 

The  program  has  one  input  card.  The  Parameter  card  specifies  the  value  for 
which  all  multiples  are  to  be  found  or  interpolated,  the  number  of  input  tape(s), 
the  number  of  files  to  be  skipped  on  the  input  tape(s),  the  type  of  data  to  be 
read  (uncorrected  fathoms  or  meters  or  corrected  fathoms  or  meters  for  bathym- 
etry or  total  or  residual  magnetic  intensity  for  magnetics),  the  specific  portion 
of  data  to  be  read,  and  the  type  of  output  desired. 

With  all  the  required  parameters  defined,  the  program  starts  to  read  the  input 
tape  one  record  at  a time.  If  only  certain  portions  of  the  data  are  to  be  used, 
the  program  will  check  to  see  that  the  fix  falls  on  the  defined  grid  or  that  it 
was  taken  on  or  between  the  two  specified  dates. 

The  program  will  continue  reading  the  first  input  tape  until  it  reads  an  end-of- 
file  mark  or  a fix  taken  after  the  last  specified  date.  If  there  are  any  additional 
tapes,  the  program  will  read  them  in  a similar  manner.  The  number  of  files  to 
be  skipped  over  on  the  new  input  tape  and  the  beginning  and  end  dates  are 
contained  on  an  Extra  card.  A maximum  of  four  input  tapes  can  be  used. 
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All  the  values,  read  or  interpolated,  will  either  be  written  on  a new  tape,  be 
printed  out  on  logical  unit  61,  or  both. 

2.2  Problem  Background 

This  program  was  written  to  aid  the  oceanographer/geophysicist  in  producing 
contour  charts  of  bathymetric  and/or  magnetics  data.  HIGHLOW  identifies 
the  high  and  low  values  and  finds  or  interpolates  values  of  a certain  interval 
(every  500  meters,  every  50  gammas).  The  program  presupposes  edited  data. 


3.0  USAGE 

3.1  Calling  Sequence  or  Operation  Procedure 
Not  applicable. 

3.2  Arguments,  Parameters,  and/or  Initial  Conditions 
Not  applicable. 

3.3  Space  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage: 

4035  octal  (2077  decimal)  locations  exclusive  of  system  library  functions. 

3.3.2  Common  Blocks: 

None. 

3.3.3  Temporary  Storage: 

None. 

3.7  Input 

There  is  one  Parameter  card  which  specifies  the  number  and  content  of  the  input 
data  tapes,  the  requirements  for  selecting  data  points,  and  the  type  of  output. 

An  Extra  card  must  be  provided  for  each  additional  input  tape.  The  actual 
bathymetric  and/or  magnetics  data  is  read  in  via  magnetic  tape  on  logical  units 
10-13.  Appendix  A presents  a sample  of  the  input  data  formats.  Appendix  B is 
a complete  description  of  the  input  deck  setup. 

3.8  Output 


The  program  will  write  all  the  required  fixes  and  values  on  the  standard  printer 
(logical  unit  61)  and/or  on  a magnetic  tape  (logical  unit  40).  The  format  of 
the  output  tape  will  be  the  same  as  that  of  the  input  tape  (Appendix  A). 
Appendix  C presents  sample  output  listings. 
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3.9  Formats 

Appendix  B describes  the  Program  deck  structure. 

3.10  External  Routines  and  Symbols 

SKIPFILE,  XABSF,  ABS,  MOD,  REWIND,  ENDFILE. 

3.11  Timing 

The  time  required  depends  on  the  number  of  data  read  and  the  number  of  fixes 
interpolated. 

3.12  Accuracy 

Not  applicable. 

3.13  Cautions  to  Users 
None. 

3.14  Program  Deck  Structure 

Appendix  B describes  the  program  deck  structure. 

3.15  Reference 

M.L.  Blodgett  and  J.V.  Massingill,  “A  Program  for  Storing  Oceanographic  Data 
on  Magnetic  Tape,”  NRL  Report  7861,  March  1975. 

4.0  METHOD  OR  ALGORITHM 


Not  applicable 


5.0  FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 
The  flow  chart  and  listing  are  given  in  Appendixes  D and  E. 

6.0  COMPARISON 


No  other  known  programs  are  available  for  comparison. 

7.0  TEST  METHOD  AND  RESULTS 

Sample  bathymetric  and  magnetics  output  listings  are  given  in  Appendix  C. 

8.0  REMARKS 


None. 
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Appendix  B 
DECK  ASSEMBLY 


7 Program  HIGHLOW  (binary  deck).  (If  the  Fortran  source  dedr  is  used  instead  of  the 
binary  deck,  a Fortran  card  is  required  after  the  Equip  card.  In  addition  a Scope  card 
and  Load  card  must  follow  the  source  deck.) 


Equip  cards  (logical  units  10. 11, 12,  13.  and  40) 


Number 

Card 

Title 

Column 

Number 

Description 

1 

Job 

1-21 

7/9  JOB,  Charge  No.,  ID  No.,  time.  See  page  2-2  of 
the  3600/3800  Computer  System  Drum  Scope 

Manual. 

2-6 

Equip 

1-18 

7/9  EQUIP,  10  = **,  RO,  HI 

7/9  EQUIP,  40  = **,  WO,  HI 

RO  = read  only; 

WO=  write  only; 

HI  = high  density. 

See  page  2-3  of  the  3600/3800  Computer  Drum 

Scope  Manual. 

7 

Program 

Deck  of 
HIGHLOW 

This  is  the  main  program.  If  the  Fortran  source  deck 
is  used  instead  of  the  binary  deck,  a Fortran  card  is 
required.  The  Fortran  card  reads  7/9  FTN,  L,  R,  X. 

In  addition  a Scope  card  with  SCOPE  starting  in 

Column  10  and  a Load  card  must  follow  the  source 
deck. 
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Number 

8 

9 


Card 

Title 


Run 


Parameter 


Column 

Number 

1-13 


1-10 


Description 

7/9  RUN,  T,  P,  R,  M,  D 
T = time  limits  in  minutes; 

P = maximum  number  of  print  or  write  operations; 

R,  M,  and  D may  be  left  blank.  See  page  2-15  of  the 
3600/380 0 Computer  System  Drum  Scope  Manual. 

1000  (IDEPTH) 

Specified  interval.  Needed  only  if  column  18  = 1. 


11  1, 2,  3,  or  4 (NOTP) 

Number  of  input  tapes.  There  can  be  a maximum  of 
four  tapes. 

12  1,2,  3,  4,  5,  or  6 (KIND) 

1 = value  in  uncorrected  fathoms; 

2 = value  in  corrected  fathoms; 

3 = value  in  uncorrected  meters; 

4 = value  in  corrected  meters; 

5 = value  is  total  magnetic  intensity; 

6 = value  is  residual  magnetic  intensity. 


14  1 or  2 (LOHI) 

1 = find  all  high  and  low  values; 

2 = do  not  find  all  high  and  low  values. 

16  1,  2,  or  3 (INEED) 

1 = write  information  on  tape; 

2 = print  out  information; 

3 = write  on  tape  and  print  out  information. 


18  1 or  2 (IMOD) 

1 = find  the  fixes  for  all  values  which  are  multiples 

of  a specified  value  (columns  1-10); 

2 = do  not  find  fixes  for  all  values  which  are  multi- 

ples of  a specified  value. 

20  2 (ISKIP) 

Number  of  files  to  be  skipped  on  the  first  input  tape. 


22  1 or  2 (IFILE) 

1 = bathymetry  file  to  be  read; 

2 = magnetics  file  to  be  read. 


7 


Blodgett 


Number 


Card 

Title 


Column 

Number  Description 


24  0, 1 or  2 (ICOPY) 

0 = only  the  data  within  a certain  area  are  to  be 

copied.  The  latitude  and  longitude  values  defin- 
ing the  area  are  in  columns  41-80. 

1 = Only  the  data  taken  during  a certain  time 

interval  are  to  be  copied.  The  beginning  and  end 
times  for  input  tape  1 are  defined  in  columns 
25-40. 

2 = All  the  data  are  to  be  copied. 

Columns  25-40  are  used  if  column  24  = 1. 

25-32  07172130 

This  is  the  time  the  program  starts  copying  the  data. 

Columns  25-26  = month; 

columns  27-28  = day  ; 

columns  29-30  = hour; 

columns  31-32  = minutes. 

33-40  08152215 

This  is  the  time  the  program  stops  copying  the  data. 

Columns  33-34  = month; 

columns  35-36  = day; 

columns  37-38  = hour; 

columns  39-40  = minutes. 

In  this  example  for  columns  25-40,  the  data  between 
July  17  and  21  hours  and  30  minutes  and  August  15 
at  22  hours  and  15  minutes  will  be  copied  from  input 
tape  1. 

Columns  41-80  are  used  if  column  24  = 0. 

41-50  40.50 

The  southernmost  latitude  of  the  area  to  be  copied. 

51-60  50.00 

The  northernmost  latitude. 

61-70  -28.00 

The  westernmost  longitude. 

71-80  -14.00 

The  easternmost  longitude. 
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Card 

Number  Title 


10  Extra 


Column 

Number  Description 

In  this  example  for  columns  41-80,  the  positions 
are  in  degrees  and  hundredths  of  a degree  (not 
degrees  and  minutes).  Southern  latitudes  and  western 
longitudes  are  preceded  by  a negative  sign  (- ). 

1-4  0, 1,  or  2 (ISKIP) 

Number  of  files  to  be  skipped  on  the  second  input 
tape. 

5-12  02250330 

Date  and  time  of  the  first  data  point  to  be  read  and 
plotted  from  the  second  input  tape. 

13-20  02280830 

Date  and  time  of  the  last  data  point  to  be  read  and 
plotted  from  the  second  input  tape. 

There  must  be  an  Extra  card  for  each  additional 
input  tape.  Since  there  is  a maximum  of  four  input 
tapes,  the  maximum  number  of  Extra  cards  is  three. 
The  dates  for  the  first  input  tape  are  on  the  parame- 
ter card. 


11  End-of-file 


Terminates  the  run. 
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Appendix  C 

SAMPLE  OUTPUT  LISTINGS 


we 

LEW 

C*T 

values 

hr 

Ml  N 

latitude 

LEkGITLDE 

depth 

IK  P6PRFCTED  *ETfcPs 

a 

IS 

10 

400 

5S.41C5 

-7.5401 

1 1 30 

H 

8 

IS 

11 

350 

59.6197 

-7.4533 

954 

L 

e 

JS 

IP 

a 

59 ,714’ 

-7,4124 

967 

H 

B 

IS 

IP 

4C0 

5S . 8654 

-7.3445 

932 

L 

8 

IS 

IP 

440 

5S . 8605 

-7,3377 

967 

H 

8 

IS 

IP 

510 

59.SC7S 

-7,3268 

789 

L 

8 

IS 

IP 

550 

59.9566 

-7,3237 

819 

H 

6 

IS 

13 

250 

60.C375 

-7,2754 

473 

L 

8 

IS 

14 

300 

60.2754 

-7,1559 

1159 

H 

6 

IS 

17 

95 

60.8445 

-6,7773 

156 

L 

8 

IS 

17 

150 

60  ,e«33 

-6 , 76 i 2 

1»3 

H 

8 

IS 

17 

2uO 

6C.8E 07 

-6,7467 

156 

l 

8 

IS 

17 

250 

60 .6981 

-6,7321 

186 

M 

8 

IS 

18 

250 

61.1157 

-6 ,5579 

130 

L 

8 

is 

19 

200 

61.3176 

-6,4789 

234 

M 

8 

IS 

19 

250 

61 .3347 

-6,4194 

221 

L 

6 

IS 

19 

350 

61.3686 

-6,4006 

336 

N 

8 

1$ 

19 

450 

61.4025 

-6,3817 

311 

L 

8 

IS 

19 

500 

61.4194 

-6,3727 

376 

N 

8 

IS 

?0 

100 

61.4914 

-6,3328 

262 

L 

8 

IS 

?0 

150 

61.5104 

-6,3725 

302 

H 

8 

IS 

PO 

POO 

61.5594 

-6,3122 

272 

L 

8 

IS 

PO 

250 

61.5485 

-6,3019 

297 

M 

8 

IS 

Pi 

575 

61 .8E7E 

-6,1208 

106 

L 

8 

IS 

?p 

350 

65.0224 

-6,0458 

1«1 

W 

8 

IS 

pp 

550 

62 .0936 

-6,0078 

97 

L 

8 

IS 

P3 

0 

62.1117 

-6,9983 

130 

N 

8 

IS 

P3 

100 

62,1503 

-5,9788 

104 

l 

8 

IS 

P3 

250 

62.1724 

-5,9675 

176 

H 

8 

IS 

P3 

550 

62.1686 

-5,9693 

101 

l 

8 

?0 

3 

50 

63.0251 

-5,7248 

1889 

H 

8 

?c 

3 

100 

63.0434 

-5,7195 

1385 

L 

8 

?c 

3 

250 

63.0986 

-5,7037 

2027 

M 

8 

70 

3 

350 

63,1356. 

-6,6910 

2016 

l 

8 

?c 

3 

550 

63.2112 

-5,6609 

P1S2 

H 

8 

?c 

4 

45 

63.2474 

-5,6497 

2108 

L 

8 

PC 

4 

ICO 

63.2687 

-5,6457 

-5,6428 

2189 

M 

8 

PC 

4 

130 

63.2603 

2167 

L 

8 

PC 

4 

200 

63.3C74 

-5,6371 

2269 

M 

6 

PC 

4 

255 

63,3586 

-5,6326 

1979 

l 

8 

PC 

4 

300 

63.346C 

-5,6789 

2001 

H 

8 

PC 

4 

450 

63,4011 

-5,6148 

-5,4757 

1970 

L 

8 

PC 

6 

550 

63.8661 

345l 

M 

8 

PC 

7 

50 

63.9222 

-5,47l3 

3137 

l 

8 

PC 

7 

360 

64.0385 

-5,4351 

3400 

M 

8 

7C 

7 

380 

64.0454 

-5,4335 

3360 

l 

8 

PC 

8 

C 

64.1517 

-5,4150 

3585 

W 

8 

PC 

8 

500 

64.3566 

-5,3677 

2590 

l 

8 

PC 

9 

0 

64.365C 

-5,3367 

3046 

H 

8 

PC 

9 

150 

64,4134 

-5,3385 

2467 

l 

8 

PC 

9 

35  0 

64.4764 

-5,3405 

3208 

H 

8 

PC 

9 

400 

64.4921 

-5,3410 

3101 

L 

8 

PC 

9 

550 

64.5394 

-5,3424 

3259 

H 

8 

PC 

10 

0 

64.5551 

-5,3429 

3188 

L 

6 

PC 

10 

50 

64 , 5709 

-5,3434 

3335 

M 

8 

PC 

10 

230 

64.6181 

-5,3449 

2493 

L 

8 

PC 

10 

230 

64.6575 

-5,3452 

2P«1 

M 

8 

PC 

10 

300 

64.6496 

-5,3459 

2544 

l 

8 

70 

10 

340 

64 .6622 

-5,3463 

3005 

H 

A 
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ALL  VALLES  WHCH  ARE 

i MUH1ELES  EF 

10 

*0 

CAY 

HR 

MIN 

LATIUCF 

LEAGIUDE 

DEPTH  IA  OERHfCTED  METERS 

B 

IS 

10 

0 

5S.2 678 

-7,5888 

105c 

8 

IS 

10 

50 

5S  .2815 

-7,5825 

1060 

e 

IS 

10 

50 

59.2421 

-7,5782 

1070 

8 

IS 

10 

ICO 

59.3C51 

-7,5736 

1080 

8 

IS 

10 

150 

59.3217 

-7,5684 

10»0 

8 

IS 

10 

150 

55.3428 

-7, 5616 

1100 

8 

is 

10 

250 

59 .3848 

-7,5546 

1110 

8 

IS 

10 

250 

59.3877 

-7,5473 

1120 

8 

is 

10 

400 

59.4105 

-7,5401 

1 1 Jo 

8 

IS 

10 

4C0 

59.4232 

-7,5360 

1120 

8 

IS 

10 

4 00 

59.4359 

-7,5320 

1110 

8 

IS 

10 

5 'JO 

59 , 4486 

-7,5279 

1100 

8 

is 

10 

500 

59 , 4!81 

-7,5249 

10»0 

8 

IS 

10 

550 

59 . 4#  7< 

-7.5219 

1080 

8 

IS 

10 

550 

59.4729 

-7,5202 

1070 

8 

is 

10 

550 

59.4782 

-7,5185 

1060 

8 

is 

10 

550 

59.4835 

-7,5168 

1050 

8 

IS 

11 

0 

59.4918 

-7,5128 

10«0 

8 

is 

11 

0 

59.5C3S 

-7,5054 

1030 

8 

is 

11 

0 

59.5142 

-7,4979 

1020 

8 

IS 

11 

100 

59.5275 

-7,4916 

1010 

8 

IS 

11 

100 

59.5368 

-7,4878 

iooo 

8 

IS 

11 

150 

59,5463 

-7,4838 

99o 

8 

is 

11 

150 

59.5589 

-7,4786 

98o 

8 

IS 

11 

230 

59.5738 

-7,4724 

97o 

8 

IS 

11 

250 

59.599C 

-7,4619 

960 

8 

IS 

11 

350 

59.6834 

-7,4344 

960 

8 

is 

12 

0 

59.7709 

-7,3870 

960 

8 

is 

12 

330 

59.8308 

-7,3601 

950 

8 

IS 

12 

300 

59.8413 

-7,3554 

940 

8 

IS 

12 

350 

59,8494 

-7,3517 

930 

8 

IS 

12 

350 

59.8551 

-7,3491 

920 

8 

IS 

12 

350 

59.8808 

-7,3466 

910 

8 

IS 

12 

400 

59.8873 

-7.34J7 

910 

8 

is 

12 

400 

59.8898 

-7,3426 

920 

8 

IS 

12 

400 

59.8719 

-7,3416 

930 

8 

is 

12 

400 

59.8742 

-7,3405 

940 

8 

IS 

12 

400 

59.8768 

-7,3395 

950 

8 

is 

12 

430 

59.8789 

-7,3384 

960 

8 

is 

12 

440 

59.8818 

-7,3373 

960 

8 

IS 

12 

440 

59.8831 

-7,3367 

950 

8 

is 

12 

440 

59.8847 

-7,3360 

940 

8 

is 

12 

440 

59.8862 

-7,3354 

930 

8 

is 

12 

440 

59.8877 

-7,3348 

920 

8 

IS 

12 

440 

59.8893 

-7,3342 

910 

8 

is 

12 

440 

59.8908 

-7,3336 

900 

8 

IS 

12 

440 

59.6924 

-7,3330 

89o 

8 

IS 

12 

440 

59.8939 

-7,3324 

8*0 

8 

is 

12 

440 

59.8454 

-7,3318 

870 

8 

IS 

12 

440 

59.8970 

-7,3311 

860 

8 

IS 

12 

440 

59,8985 

-7,3305 

850 

8 

is 

12 

440 

59.900C 

-7.3299 

840 

8 

is 

12 

440 

59. 9014 

-7,3293 

830 

8 

is 

12 

440 

59.9031 

-7,3287 

820 

8 

1$ 

12 

440 

59.9047 

-7,3281 

810 

8 

IS 

12 

440 

59.9062 

-7,3275 

800 

8 

is 

12 

440 

59.9077 

-7,3269 

790 

8 

is 

12 

510 

59.9085 

-7,3267 

790 

8 

is 

12 

510 

59.9148 

-7,3257 

BOO 
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h!UH  »NC 

iew 

VALUES 

and  all 

VALLES  kHIO 

ARE  HLLTIHLES 

10 

*0 

CAY 

HR 

MIN 

LATIUCF 

lekoitlde 

DFPTh  D rCRHtCTEI  “ETfcNS 

a 

IS 

10 

s 

5S.2f7e 

-7,5888 

1050 

e 

IS 

10 

50 

59,2815 

-7,5625 

1 0*0 

8 

IS 

10 

50 

59.2921 

-7,5782 

1070 

a 

IS 

10 

100 

59 . 3C 5 1 

-7,5736 

1080 

a 

IS 

10 

150 

59 . 3? 1 7 

-7,5484 

1 0»0 

a 

IS 

10 

150 

55,3426 

-7,5416 

1100 

a 

IS 

10 

250 

59.3848 

-7,5546 

1110 

a 

IS 

10 

250 

59.3877 

-7,5473 

1170 

a 

ts 

10 

4-0 

59.4105 

-7 ,5401 

1130  H 

a 

IS 

10 

400 

59.4732 

-7,5360 

1120 

a 

IS 

10 

400 

59.4J5S 

-7,5320 

1 110 

a 

IS 

10 

5G0 

59.4488 

-7,5279 

1100 

a 

IS 

10 

500 

59.4581 

-7,5249 

1090 

a 

IS 

10 

550 

59 .48  78 

-7,5219 

10*0 

a 

IS 

10 

550 

59.4729 

-7,5202 

1070 

a 

IS 

10 

550 

59,4782 

-7,5185 

10*0 

a 

IS 

10 

550 

59.4835 

-7 ,5160 

1050 

a 

IS 

11 

0 

59.4918 

-7,5128 

10*0 

a 

IS 

11 

0 

59 .5035 

-7,5054 

10  50 

a 

IS 

11 

0 

59.5162 

-7,4579 

1020 

8 

IS 

11 

1 JO 

59.5775 

-7,4516 

1010 

a 

IS 

11 

130 

59,5368 

-7,4678 

1000 

a 

IS 

11 

150 

59.5463 

-7,4838 

9*0 

a 

IS 

11 

150 

59.5585 

-7,4786 

9*0 

a 

IS 

11 

230 

59.5738 

-7.4724 

97o 

a 

IS 

11 

250 

59.5590 

-7,4619 

9*0 

a 

IS 

11 

350 

59.6197 

-7,4533 

954  L 

a 

IS 

11 

350 

59.6834 

-7,4344 

9*0 

8 

IS 

12 

0 

59.7143 

-7,4124 

9*7  H 

a 

IS 

12 

0 

59.7705 

-7,3670 

9*0 

a 

IS 

12 

300 

59.8308 

-7,3601 

950 

a 

IS 

12 

300 

59,8413 

-7,3554 

940 

a 

IS 

12 

350 

59.8494 

-7,3517 

930 

a 

IS 

12 

350 

59.8551 

-7,3491 

920 

a 

IS 

12 

350 

59.8808 

-7 , 3466 

910 

a 

IS 

12 

400 

59.8854 

-7,3445 

902  l 

a 

IS 

12 

400 

59.8673 

-7,3437 

910 

8 

IS 

12 

400 

59.8694 

-7,3426 

920 

a 

IS 

12 

400 

59.8719 

-7,3416 

930 

a 

IS 

12 

400 

59.8747 

-7,3405 

940 

a 

IS 

12 

400 

59.8766 

-7,3395 

950 

a 

IS 

12 

400 

59.8785 

-7,3384 

9*0 

a 

IS 

12 

440 

59.8805 

-7,3377 

9*7  H 

a 

IS 

12 

440 

59.8616 

-7,3373 

960 

a 

IS 

12 

440 

59.8831 

-7,3367 

950 

a 

IS 

12 

440 

59.8847 

-7,3360 

940 

a 

IS 

12 

440 

59.8862 

-7,3354 

930 

a 

1« 

12 

*40 

59.8677 

-7,3348 

920 

a 

IS 

12 

440 

59.8693 

-7,3342 

910 

a 

IS 

12 

440 

59,6908 

-7,3336 

900 

a 

IS 

12 

440 

59.8524 

-7,3330 

890 

a 

IS 

12 

440 

59.8539 

-7,3324 

880 

a 

IS 

12 

440 

59.8554 

-7,3318 

870 

a 

IS 

12 

440 

59.8570 

-7,3311 

6*0 

a 

IS 

12 

440 

59.8585 

-7,3305 

650 

a 

IS 

12 

440 

55.9000 

-7,3299 

840 

a 

IS 

12 

440 

59.9014 

-7,3293 

B30 

a 

IS 

12 

440 

59,9031 

-7,3287 

820 

a 

IS 

12 

4*0 

59.9047 

-7,3281 

BIO 
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Appendix  E 

SOURCE  LANGUAGE  LISTING 


PROGRAM  HIGHLOW 

REAL  LATMIN.LATMAX.LONGMIH.LONGMAX 
0 IMENS I ON  XLATC3),XLONGC3>,IREADC3>,READ(2> 

DIMENSION  M0NTHC3),IDAYC3).HRC3),XNINC3> 
dimension  MtSSAGEC2> 

0ATACMESSAGE=2H  L,2H  H) 

ILBL=C 
I EX=0 
IT  = l" 

IFUOGE  = G 

C READ  PARAMETER  CARD 

lOGO  PEADC6C ,25O)IUEPTH,N0TP,K INO.LOHI, INEED,  IM  10.1SKIP 

I,IFILE.IC0PY,IDATE1,IHR1,IDATE2. IMR2 , L AT  N I N.  L AT  M AX  .LONGM  l N,  L0NGM  AX 

2 5G  FORMAT  cm  ,211 , 612, 515, 5F 10. 5) 

DEPTH= IDEPTH 
5C7  IREC=l 
NUM  = 1 

DO  335  1 = 1 , I SKIP 
CALL  SKIPFILFCIT) 

3 34  CONTINUE 

IFCIFILE.E3.2)  GO  TO  60u 

C READ  8ATHYMCTRY  DATA 

20k  RE AUC IT, 38)1  YEAR ,MONTH( 1),IDAY(1),HR(1),XMIN(1),XLAT(1 ), XL6NGC I ) • 
1UNF ATM, ICKMET 

38  FORMAT < 1 3X312. 1XF2.F3.F8.5.F9.5, IJXF5, I5.13X) 

IFC IOCMECK, IT ) 2C  o , 55 
55  IFCfcOF,IT)SjO,SU 


C CHECK  IF  OATA  TO  BE  COPIED 

50  UNF  S = UNF  ATH 
ICP  S= 1C  RME  T 

IFC  IC OPT- 1)02 5.926.92 7 

925  IFCXLATCD.LT.LATMIN.OR.XLATCD.GT.LATMAX.  JR . XL ONGC  I ) . L T . LONGM  IN 
l.OP.XLPNGCD.GT.LONGMAX)  GO  TO  20'. 

GO  TO  927 

926  IDAYI=»"NTHC 1 )*10C  ♦ IDAYC1) 

MINUTE=HRC1 )*10C  ♦ XMINC 1 )/ 10 
IFCIOAYl.LT. IUATEI)  GO  TO  2uP 
1FCIDAY1 ,GT. IOATE2)  GO  TO  5G0 

IFCIOAY1.EQ.  IOATE1.  ANU.MINUTr.LT.  IHRDGO  TJ  20" 

IFCIOAVl.t'J.  IU ATE 2.  AND. "I NUT E .GT  . IHR 2 )GO  T1  5Crt 

C CONVERT  DEPTH  DATA  IF  REQUIRCD 


IBIS  PAO*  IS  BSSI  QVULin  FRAOIiOAW** 
fjyjH  Q&n  WHBISH®  TO  ECO 


14 


NRL  REPORT  8305 


927  IFCKIND.FQ.5)  00  TO  55 
IFCKINf)-2)l01»ll«2«ll3 

1P1  IREA0C1  )=UNFATll/10.U 
GO  TO  105 

lr  2 READC1 >=ICRMET/1 .8288 

IF ( RE AU(l). LT. 0)00  TO  551 
IREADC 1 )=  RtAUCl  >*.5 
GO  TO  105 

551  IREAOC  1 >=  REA0(l)--.5 
GO  TO  105 

10  3 RE  ADO  > = (UNFATH/10.")»  1.8  2 88 

IF ( REA0C1 ).LT.O)GO  TO  552 
IREAOC 11=  READCl )*.5 
GO  TO  1 C 5 

552  IREADC 1 ) = RtA0Cl)-.5 
GO  TO  105 

C READ  MAGNE  TICS 

6<K>  READCIT,6.5)IYEAR,M0NTHCl),IUAYCl),HRCl>,X  UNCI  > , XL  AT C 1 ) , XL ON& C 1 > , 
1 ITM I , I RM I 

6r*  5 FORMATE  13X312,  lXF2,F3,r  8.5,  F9.5, 23X2  15,5  X) 

IFCIOCHECX,IT)60O,6Cl 

60 1 IFCEOF, IT)500,602 

C CHECK  IF  DATA  TO  PE  COPIED 

602  ITM IS= I TM I 
IRMIS=IRMl 

IFCICOPY-1 >928,929,930 

928  IFCXLATCD.LT  .LATMIN.OR.XLAT  Cl  ).GT  .LATMAX.  JR  . XL  ONGC  1 > . L T . LON&M  IN 
1 .OP.XLONGC 1 J.GT.LON&MAX)  GO  TO  600 

GO  TO  930 

929  IDAYDMPNTMC1  >»1C0  ♦ IDAYCl) 

KINUTE=HRC1>*100  ♦ XMINC1)  / 10 
IFCIDAY1.LT. IDATEl)  GO  TO  60" 

IFC  IOAY1 .GT. IOATE2)  GO  TO  500 

IFCIDAY1.E0.TUATEI.AN0. MINUTE. LT.IHRDGO  TJ  6C>9 
IFCIOAY1 ,E0. IOATE2. AND. MINUTE .GT . IHR2)GO  TJ  509 
9 31'  . IFCKIND.E0.5)&0  TO  603 
IRE ADC  1 )=IRMI 
GO  TO  105 

603  IREAOC 1 ) = I TM  I 
GO  TO  1C  5 

55  IREADC 1 )=ICRMET 

C CHECK  IF  MULTIPLES  OF  SPECIFIED  DEPTH  REQUIRE) 

105  IFCIMOD-l)5u2, 305,502 
C WRITE  LABELS  AT  TOP  OF  PAGE 


BLODGETT 


888  IFELOHI.EQ.l . AND.  IMOO.E'J.l  ) CO  T9  969 
IFELCHI.  EQ.l.  AND.  IMOO.  EQ.2>  GO  TO  910 
WRITEE61.913)  IDEPTH 

913  FORMATE 1H1 ,33HALL  VALUES  WHICH  ARE  MULTIPLES  OF, IIP) 

GO  TO  914 

909  WRITEE61.911)  IDEPTH 

911  FORMATC1H1 .57HHIGH  AND  LOW  VALUES  ANO  ALL  VALUES  WHICH  ARE  HULTIPL 
1ES  OF, IIP) 

GO  TO  914 

910  WRITEE61 ,912) 

912  F0RMATC1H1.19HHIGH  AND  LOW  VALUES) 

GO  TO  914 

914  GO  TOE901, 962, 903, 90 4, 915,916), KIND 
9P1  WR I TE ( 6 1 , 90 5 ) 

905  FORMATE 1H  ,88H  MO  DAY  HR  MIN  LATITUOE  LONGITUOE 

1 DEPTH  IN  UNCORRECTEO  FATHOMS/) 

GO  TO  164 

9P2  WRITEE61.906) 

906  FORMAT E 1 H ,86H  HO  DAY  HR  MIN  LATITUDE  LONGITUOE 

1 DEPTH  IN  CORRECTED  FATHOMS/) 

GO  TO  164 

9n  3 WRITEE61.90T) 

9P7  FORMATE  1 H ,87H  MO  DAY  HR  MIN  LATITUDE  LONGITUDE 

1 DEPTH  IN  UNCORRECTEO  METERS/) 

GO  TO  164 

915  WRITEE61 ,91 7) 

917  FORMATE  1 H ,84H  MO  DAY  HR  MIN  LATITUDE  LONGITUDE 

1 TOTAL  MAGNETIC  INTENSITY/) 

GO  TO  164 

916  WRITEE61 ,918) 

918  FORHATC1H  ,87H  MO  DAY  HR  MIN  LATITUOE  LONGITUDE 

1 RESIDUAL  MAGNETIC  INTENSITY/) 

GO  TO  164 

904  WR ITE  C 6 1 > 9P  8 ) 

9P8  FORMAT E 1H  ,85H  MO  DAY  HR  MIN  LATITUDE  LONGITUDE 

l DEPTH  IN  CORRECTED  METERS/) 

164  NUM-1 

IFEITRIP.EQ.4)  GO  TO  472 
IFEITRIP-2)  162,172,1602 

162  WRITEE6l,4r)M0NTHEl),IDAYEl),HREl),XMINEl),XLATEl)«XL0NGEl),IREAD 
1 E 1 ) 

4C  FORMATE 1H  ,5XI2,I4,4XF2,3XF3,5XF8.4,5XF9.4,10XI5,1XA2) 

NUM  =NUM  ♦ 1 
GO  TO  3C5 

305  IFEIFILE.EQ.2)  GO  TO  605 
C READ  BATHYMETRY  DATA 

315  REAOEIT,38)IYEAR,MONTHE2),IDAYE2),HRE2),XMINE2),XLATE2),XL0NGE2), 
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1UNFATH, ICRPE T 
IF( IOCHECK  * IT)  315, 70u 
70t  IFCE8F«IT35uO,701 
7n 1 IFCICOPY-13931, 932, 3333 

C CHECK  IF  DATA  TO  8E  COPIED 

9 31  IFCXLATC23.LT.LATPIN.OR.XLATC23.GT.LATPAX.  7R  .XL0NGC2 3 . LT  . LONGP IN 
1 .0R.XL0NG(2).GT.L0NGHAX)  GO  TO  315 
GO  TO  3333 

932  IDAV1=P0NTHC2 )» IOC  ♦ IDAYC2)  1 

P INUTE  = HRC  2 3 * IOC  ♦ XPINC23/10 
IFCI0AV1.LT. IDATEi)  GO  TO  315 
IFCIDAV1.GT.IDATE2)  GO  TO  500 

IFCIDAY1.EQ.I0ATE1. AND. PINUTE.LT. I HR 1)G0  TT  315 
IFCIOAY1.EQ.IOATE2.ANO.PINUTE.GT.IHR23GO  T1  50? 

GO  TO  3333 

C CONVERT  DEPTH  DATA  IF  REQUIRED 

3333  IFCKINO.EU.63  GO  TO  705 
IFCKIN0-237Q2, 703,706 

702  IREADC23=UNFATH/10.0 
GO  TO  756 

703  RE AO C23=ICR MET/1.82 88 

IFC  READC23.LT.ilGO  TO  553 
IREADC  2 3=  REA0C23  + .5 
GO  TO  756 

553  IREADC2  3=  READC23-.5 
GO  TO  756 

7f 6 RE  ADC 2 3 = C UNF  ATH/10 , 03*1.8288 

IFC  READC23.LT.u3GO  TO  556 
IREADC  2 3=  READC23+.5 
GO  TO  756 

556  IRE AD C2 3=  REAOC23-.5 
756  IUET=XA8SFC1READC13-IREADC233 
GO  TO  706 

C READ  PAGNETICS 

605  READCIT,6063IYEAR,PONTHC23,IDAYC23,HRC23,XMINC23,XLATC23,XLONGC23, 

1 ITPI , IRPI 

IFC IOCHECK , I T 360  5 ,6«6 

606  IFCEOF,IT15uO»607 

C CHECK  IF  DATA  TO  BE  COPIEO 

607  IFCICOPY-13 9 36,935, 936 

936  IFCXLATC23.L7.LATPIN.OR.XLATC23.GT.LATPAX. JR . XLONGC2 3 . L T . LONGP IN 
1.0R.XL0NGC23.GT.L ONGPAX 3 GO  TO  605 


17 


BLODGETT 


GO  TO  936 

935  IDAY1=MONTH<23*100  ♦ IDAYC23 
MINUTE =HR{ 23* 100  ♦ XMIN<2)/10 
IFTIOAYl.LT. I0ATE1)  GO  TO  605 
IFUDAY1.GT.IDATE2)  GO  TO  500 

IFCIDAYl.EQ.IOATEl. AND. MINUTE. LT.IHR13GO  TJ  6C5 
IFU0AY1.E0. IOATE2. AND. MINUTE. GT. IHR23GO  TJ  500 

936  IFCKIN0.EQ.5)  GO  TO  608 
IREAOT  2 ) = IRNI 

GO  TO  746 

608  IREADC23-ITMI 
GO  TO  746 

705  IREADC2 3=ICRMET 
GO  TO  756 

746  IU8T=XA8SF(1READ(1 3-IREADC2 33 
C CHECK  IF  ANY  MULTIPLES  OF  SPECIFIED  DEPTH  BETTE EN  TWO  FIXES 

706  SLCT1=IREAD<1 3/  DEPTH 
I SLOT  1 = SLOT  1 
SLOT2=IREAOC23/  OEPTH 

1 SL  OT2  = SLOT  2 

IFCHOD(IREAD<13»IOEPTH3.EQ.l3GO  TO  97C 
ILHL=0 

IFC ISL OT 2- I SLOT1 3900,833,301 
444  XLATTl 3=XLAT<23 

XLONGC1 3=XLONG<23 
IRE  AO( 1 J = IREAD<23 
MONTH C 13=MONTH<23 
IDAYC1 3-IDAYC23 
HRC13=HRC23 
XMINC1 3-XMINC23 
GO  TO  305 

90C  IFCIREAOd 3.LT.C . AND. I RE AO< 2 3 .LT . 0 3 GO  TO  935 

2 NUM  = I SLOT1-SLOTI 
I SLOT  1 = I SL  OT  1 - 1 
GO  TO  301 

3Cu  IF(IREAD(13.GE.0.AND.IREADC23.GE.O3  GO  TJ  J„9 
ZNUM=ISLOT 1-SLOT1 
I SL  OT  1 = I SLOT  1 ♦ 1 
GO  TO  301 

905  ISLOT1 = I SL  OT 1-1 

ZNUM=ISLOT1-SLOT1 
GO  TO  301 

309  I SLOT1  = I SL  OT  1 ♦ 1 
ZNUM= I SL  OT 1 - SL  OT 1 

C INTERPOLATE  FIXES 

301  DENON=  SLOT2-SL  OT  1 


NRL  REPORT  8305 


CALAT=XLATCl 3 *C ZNUH/DE N0M3*C XL AT C 2 3-XL AT C l 33 

CAL0NG=XL0NGC1 3* CZNUM/DENOM 3*C XLONGC 2 3-XLO NGC 1 3 3 

TEST=2NUM/DEN0N 

XYZ=TEST»CIREADC23-IREADC133 

XX=IRE A0C13*XYZ 

IKXX.LT.c)  GO  TO  559 

XX=XX  ♦ .5 

GO  TO  563 

559  XX=XX  - .5 

563  IC ALHE  T = XX 

IFCICALMET.EQ.IREADCI33GO  TO  307 
IFCICALMcT.EQ. IREADC233GO  TO  307 
C PRINT  AND/OR  WRITE  ON  TAPE  THE  REQUIRED  INFORMATION 

323  IFCINEEO-2370,72,71 

7C  CO  T0C170, 170, 273, 270, 3370,9703, KIND 

17C  WR I TEC  90.171  3 I YE AR , HONTHC 1 ),  IDAYC1  3.HRC1  3.XMINC1  ) ,CAL  AT  . C ALONG . ICA 
1LNET 

171  FORMAT ( 13X ,3I2,lXF2,F3,F8.9,F9.9,lOXI5,23X3 
GO  TO  307 

270  WRITEC90,271)IYE»R.NONTHC1 3 , IDA Y C 1 3 , HRC 1 3 , XM INC l 3 ,C AL A T ,C ALONG  , IC ALHE T 
1LNET 

271  FORMATC13X3I2,lXF2,F3,F8.9,F9.9, 15X15. 18 X) 

GO  TO  30  7 

3 370  WRITEC90.3371 )I YE AR , HONTHC 1 ) . ID  AY  C 1 3 ,HR  C l 3,  XHI N C 1 3 ,C  AL  A T ,C  AL  ONG, 

1 ICALMET 

3 371  FORNATC  13X312 , 1 XF2,«=3  , F3.  A,  F9. 9 , 2 3X15 , 13X  3 
GO  TO  307 

9TL  WRI TEC  90,9 71 3 1 YEAR, HONTHC 1 ) , IDA Y C 1 ) , HR C 1 ) , XN INC  1 3 , C AL AT ,C ALONG , IC A 
ILMET 

971  FOPHATC13X3I2.1XF2,  F 3 . Fe . 9 ,F9 . 9 , 28X 1 5 ,5 X 3 
GO  TO  30  7 

71  GO  TflC57C,570, 670,670,770,8703. KIND 

77l  WRITEC90.7  71 3 I YE AR .HONTHC 1 3 , IDAYC1 3 , HRC 1 3 , XM INC  1 3 ,C AL AT ,C ALONG , IC AL 
1LHET 

771  FORHATC13X3I2,1XF2,F3,F8.9,F9.9,23XI5,10X3 

GO  TO  72 

5 70  WRITEC90.5  71  3 1 YE  AR  , HONTHC  1 3 , ID  AY  C 1 3 . HR  Cl  3 , XM  INC  1 3 ,C  AL  A T , C ALONG . IC  ALHE  T 

ILMET 

571  FORMATC13X3I2,lXF2,F3,F8.9,F9.9,lUXI5,23X3 

GO  TO  72 

6 7(i  WRITEC  9C  ,6  71  3 I YEAR  .HONTHC  1 3,  IDAYCl  3,  HRC  1 3,  XXI NC  13  , C AL  A T , C ALONG  » IC  A 

ILMET 

6 71  FORMAT C 13X312, 1XF2.1XF 2, F3.F8, 9, F9 . 9 , 1 5X  1 5 , 1 8X3 
GO  TO  72 

870  WRITEC90,8713IYEAR,HONTHC13,IOAYC13,HRC13,XMINC13,CALAT,CALONG,ICAL 
ILMET 

871  F ORMATC13X, 312, 1XF2.F3.F8. 9, F9. 9, 28X15,5X3 
GO  TO  72 

72  IFC  MODCNUM.603.NE.13  GO  TO  172 
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ITR  IP  = 2 
GO  TO  8RB 

172  WRITEC61  ,69)MONTH(l>,IDAY(l  ) , HR  ( 1 ) ,X  M I N(  l ) ,C  AL  A T ,C  ALO  NG  , I C ALN  E T 
NUM=NUM  ♦ 1 
GO  TO  307 

307  IFCIFUDGE.EU.  1 ) GO  TO  6666 

1F( ISLOT2-ISLOT1  ) 900,633,305 

6C7  IREC-2 

GO  TO  6U1 

813  IF(MOD(IREAO(2),ID£PTM>.EC,|.c)  GO  TO  333 

IFCIREADCl  >.LT.'j.ANO.  I p E ADC  2 ) . G T . C ) GO  TO  150 
IFCIREADCl). GT. 6. AND. IREAO(2).LT.G)  GO  TO  169 
GO  TO  333 

15C  IFUDGE -1 
GO  TO  3»,0 

16C  IFUDGE  = 1 
GO  TO  900 

5665  IFUDGE -0 
GO  TO  333 

975  IFCILBL.EQ.l)  GO  TO  1 
IF(INEED-2)37u.  372,371 

1 ILBL=5 

GO  TO  307 

370  GO  TOC137J, 1370, 1371, 1371, 1372, 1373), KINT 

1 37C  WRITE  0 60, 1 3 76 > I YE AR , MONTHC 1 ), IDAYCl),HRCl  I.XMINC 1 >,XLA1(  D.XLONG 
1 Cl). IRE  AO (1) 

1376  FORMAT C 13X312, 1XF2.F 3, F 9.6 ,F9. 6, 10X15,2 3X) 

GO  TO  3E  7 

1371  WRITE(60.1375)IYEAR,MONTH(l),IOAY(l),HR(l),XMIN(l),XLAT(l),XLONG(l) 
I)iIREAD(I) 

1 375  FORMAT ( 1 3X 312, 1XF2,F3,F8. 6, F9. 6,15X15,18X3 
GO  TO  3C7 

1 372  WRITE (60 ,1  376) I YEAR .MONTHC 1 ) , ID A Y ( 1 ) , HR ( 1 ) ,X M INC  1 ) ,XL A T < 1 ),XLONG( 1 ) . 

I ) ,IREAD(1) 

1376  FOPM4T(13X3I2,lXF2,F3,Ffc.6,F9.6,23XI5,lCX) 

GO  TO  3(  7 

1 37  3 WRITE (60, 1 377)1  YEAR, MONTHC 1), IUAYC1 ), HR  Cl  ), XM INC  1 ) , XL  AT C 1 ) , XL ONGC 1 ) 

II) ,  IRCAUO  ) 

137  7 FORMAT! 1 3X  312, 1XF2  »F3,FH.6 , F 9. 6 ,26X15 ,5X  ) 

GO  TO  3G  7 

371  GO  TOC2371,2371,23T2,237?,2373,  2376), CIND 

2 371  WR I TEC  6i  ,2  375  )I  YEAR,  MONTHC  1 ),  IDA  Y(1  ),HR(1),XMIN(1  ),XLAT(1  ), XL  ONoCD.IREA 

1 ) , IRE  AD  C 1 ) 

2375  FORMAT(13X3I2,lXF2,F3,F8.6,F9.6,luXI5,23X) 

GO  TO  372 

2 37  2 WRITE  (61,  ,2  376  ) IY  E AR  .MONTHC  1 ) , IDAY(1),HR(  l),XMINC  1),XLAT(  1), 

1XLONGC1 ) , I Rc  AD( 1 ) 

2 376  FORMAT!  1 3X  3 1 2 , 1 XF2  ,F  3 , F 8 . 6 , F9 . 6 , 1 5 X 15 , 1 8 X ) 

GO  TO  372 


tHIS  PAUE  IS  BSSI  ftUALin 

fhOM  Oori  /UWU.SH®  TO  ' 
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2 37  3 WRITE (*0 ,2377) I YE AR , MONTHC 1 ) , ID flY ( 1) , HR ( 1 ) , XNIN ( 1 3 , XL  ft T ( 1 ) , XL 3NGC 1 ) 
l)tlREAO(l) 

2 37  7 F0R*ATC13X3I2,1XF2,F3.F8.4,F9.4,2  3XI5,DX ) 

CO  TO  372 

2374  WRITEC40i2378)IYEAR>MONTHCl)tIDAYCl>iHRCl)tXMINCDtXLATCD«XLONCCl)» 
1 )»IREAD(1 ) 

2378  FOKHATC13X3I2,lXF2,F3tF8.4tF9.4.29XI5,5X> 

GO  TO  37 2 

372  IFCM30CNUN,60 >.NE  .1)  GO  TO  472 
I T P.  I P = 4 
GO  TO  888 

472  WRITE(61,4C)HONTHO),IDAYC1),HR(1),XMIN(I),XLATC1),XLONG<1),IREAU(1) 
11  > 

NUN=NUH  ♦ 1 
GO  TO  3C7 

C CHECK  IF  HIGH/LOW  VALUES  WANTED 
333  IFCL3HI.NE.DG0  TO  444 

401  XLONGC  3 3 -XLONG C 1 ) 

XLATC33=XLATC13 
IREADC33MRLAOC13 
MCNTHC3  3=«0NTHC1  ) 

IUAYC33=IDAYC13 
HRC  33  = HR  C 1 3 
XNINC33-XNINC1 ) 

XLATCl 3 - XL  AT (2) 

XLONGC 1 )=XLONGC2) 

IRE AD (1 3=IREADC23 
MOtUHCl 3=MONTHC2) 

I0AYC1 3=IDAYC23 
HR  Cl 3 = HR  C2  > 

XH INC  1 3 = XH INC7  3 

402  IFCIFILE.E0.2)  GO  TO  6IC 

C READ  BATHYMETRY  DATA 

412  R£A0CIT,38)lYEAR,M0NTHC2)t ID A Y (2 3 . HR C 2 ) . XMINC 2 3 , XL ATC 2 ) , XLONGC 2 ) , 
1UNFATH, ICRMET 
IFCIOCHECK,lT3412,4l'3 

403  IFCEOF,IT)5>/0,4L4 

C CHECK  IF  DATA  TO  Be  COPIEO 

404  IFCIC0PY-1)937,938,939 

937  IFCXLATC23.LT.LATMIN.OR.XLATC2 3.  GT  .L  A THA X . OR  . XL  0NGC2  3 . L T . LONGM  IN 
1 .OR.XLONGC2).GT.LONGMAX)  GO  TO  412 

GO  TO  939 

938  IDAYI=MONTHC2)*10l  ♦ IDAYC2) 

MINUTE=HRC2)*10&  ♦ XMINC23/10 
IFCI0AY1.LT.  IOATED  GO  TO  412 
IFCI0AY1.GT.  I D AT E 2 ) GO  TO  5C0 
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IFCIDAYI.EQ.IDATEl. AND. MINUTE. LT.IHR13G0  TT  612 
IFCI0AY1.EQ. IDATE2. AND. MINUTE. GT.IHR23GO  TJ  5C9 
GO  TO  939 

C CONVERT  OEPTH  DATA  IF  REOUIPED 

939  IF  CK IND. EQ . 6 3 GO  TO  805 
IFCKINO-23  802.803.806 

802  IREADC23=UNFATH/10.0 
GO  TO  826 

803  READC23=ICRMET/1.8286 

IF ( REAOC23.LT.03GO  TO  855 
IRE ADC 2 3 = R£ADC23*.5 
GO  TO  826 

555  IRE AOG  2 3=  REAO(2)-.5 
GO  TO  826 

80  8 KEAOC2 3=CUNFATH/10. 0*1.8288 
IF(  REAO(2).LT.O)GO  TO  556 
IRE ADC  2 3 = READC2 3*. 5 
GO  TO  826 

556  IRE  ADC  2 3 = REAOC23-.5 

826  IU6T=XA8$FCIREAOC13-IREAOC233 
GO  TO  806 

C READ  MAGNETICS 

6 1C  PEADCIT,6363IYEAR,MONTHC23,IOAYC23.HRC23.X1INC?3.XLATC23,XLONGC23.ITNI,IRN 
1 I TM I.  IRMI 

IFCIOCHECK.IT3610.611 

611  IFCE0F.IT3  50 0 .612 

C CHECK  IF  DATA  TO  BE  COPIED 

612  IFCIC0PV-1396J.981.962 

96C  IF  CXLAT C2  3.LT. CATKIN.  OR. XLATC2 3.  GT.LATPAX.  1R  . XL ONGC 2 3 . L T . L ONGM IN 
1.0R.XL0NGC23.GT . LONGMAX  3 GO  TO  61c 
GO  TO  962 

961  IUAY 1 = MONTHC 2 3 *100  ♦ IUAYC23 
MINUTE=HRC23*10C  ♦ XM1NC23/10 
IFCIDAY1.LT.IUATE13  GO  TO  6lt 
IFCI0AY1.GT. IOATE23  GO  TO  5G0 

IFCIDAY1 .EQ. IDATE 1. AND.MINUTE.LT. IHR1 3G0  T3  61C 
IFCIDAY1 .EO. IUATE2. AND. MINUTE. GT. IHR23G8  TJ  5CJ 
GO  TO  96? 

962  IFCKIND.EQ.53  GO  TO  613 
IREADC23=IRMI 

GO  TO  836 

613  IRE  ADC  2 3 = I TM  I 
GO  TO  "36 

805  IREA0C2 J=ICRM£ T 
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GO  TO  326 

8 36  IUBT=XAUSFC  I R E AD  C 1 ) - I K E AO  C 2 ) ) 

C CHECK  FOR  LOW  AnO  HIGH  VALUED 

8f6  IF( IRE ADC1 ).GT. I R E AO  C 3 ) . ANO  . IREAOE 1 ).GT.  IKEA0C2) )G0  TO  82G 
ACS  IFCIREADCl).LT.IREADC3).ANP.IREA0Cl).LT.  IREA0C2))  GO  TO  830. 

812  1FC  IMOD-  1 ) Ac  1 f 70  6,  A j 1 

C PRIM  AND/Ok  WRITE  ON  TAPE  THE  RE3UIKE0  INFOk CATION 
811  IF(INEE0-2)668,662,661 

668  GO  TOC  1668 , 1660 . 1661 . 1661 , 1662 , 1663). K INU 

1660  WRITE ( AO. 166ATI YEAR .MONTH ( lT.IOAYClT.HRCD.XMI N( 1 ) ,X L AT ( 1 ) , X L ONGC 
11T.IRE AGC1) 

166A  FCRMATC13X3I2.1XF2,F3,FR.A,F9.A,i0XI5,23X) 

GO  TO  812 

1661  WRITE  C At , 1665  >1  YE Ak.MONTHC 1 ) , 10  AY ( 1 ) , HR ( 1 ), KM I N ( 1 ) , X L AT ( 1 ), XL ONGC 
11). IRE  AD  Cl) 

1665  FORMATC13X3I2, 1XF2.F3.FR. A. F9. A, 1SXIS.13X) 

GO  TO  812 

1662  WRITE CAO , 1 666 ) I Y E AR , MONTHC 1 ) . 10  A Y C 1 ) . HR C 1 ) . XM I N C 1 ) , X L A T C 1 ) . XL ONGC 

1 1 ). IRE  ADC  1 ) 

1666  FORMA  TC  13X3I2,1XF2,F3,F8.A,F9.A,23XI5,  MX  ) 

GO  TO  812 

1663  WKITECAO,1667)IYEAR,MONTHC1),IDAYC1),HRC1).XMINC  1 ) , X L A T C 1 ) ,X L ONG 
lClT.IREAOCl  ) 

1667  FORMA TC 13X3 12.1 XF2.F3.F8, A .F9.A ,28X15, SX) 

GO  TO  812 

661  GO  T0C2663, 26 60,1661, 2661, 2662, 2663), KINO 

2660  WRITECA  C , 2 66  A ) I Y E AR  , MONTHC  1 ) . ID  A Y C l ) . HR  C 1 ) . X N IN  C 1 ) , XL  A TC  I ) .XL  ONGC  1 ) , I 
1), IRE ADC  1) 

266A  F0PMATC13X3I2,1XF2,F3,F8.A,F9.A, UXI5.23X) 

GO  TO  662 

2661  WR I TEC  AC;  ,2665  ) I YEAR  , MONTHC  1 ) , 1 0 A Y C 1 ) , HR  C 1 ) , X M INC  1 ) . XL  AT  C 1 ) • XL  ONGC  1 ) , 

1 T.IREAOCl  ) 

2 665  FORMAT  Cl  3X312, 1 XF2 . F3 , F8 . A , F9 . A , 15X15, 18X) 

GO  TO  66? 

266  2 WRITECAl ,2 666 ) I Y EAR , MONTHC 1 ) , I DA Y C 1 ) . HRC 1 ) .XM1NC 1 ) . XL  A T C 1 ),XL0NGC1 ), 
1), IRE  AD  Cl) 

2 666  F0RMATC1 3X 3 £2 , 1 XF2 ,F3 ,F8 . A , F9 . A , 2 3X 15 , 10X) 

GO  TO  662 

266  3 WR I TEC  AO .2667 ) I YE AR , MONTHC 1 ). 10 AY  Cl ),HRC l ).XMINC 1 ).XL ATC 1 ). XL ONGC 1 ) 
D.IREA0C1) 

2 667  TOP. MATC13X3I2,1XF2»F3,F8.A,F9.A,28XI5,5X) 

GO  TO  662 

662  IFCMOOCNUM,60).Nt.l>  GO  TO  1602 
ITR IP  = 3 

GO  TO  888 

1 6c 2 WKITEC6 1 ,A0)MONTHC 1), 1DAYC 1 ) , HR C 1 ) , X M I NC 1 ) .XL AT C 1 ) , XL ONG C 1 ) , I R E ADC  1 ) , 
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1 1 3«MrSsAGtELLL  ) 

NUM^NUM  ♦ 1 

GO  T8  £12 
b2v  LLL  -2 
I L6L=1 
GO  TO  hi* 
b3l  LLL  = 1 
I LLL= 1 
GO  TO  “lr 
5CC  I E X = 7 c X ♦ 1 

HEWING  IT 

IFENOTP .1 y. IFX 3 GO  TO  Sw 5 
I T = IT  «■  1 

RLAUE6'  ,5C6  ) I SKIP  , ID ATE  1 , IHR1  , I U ATE  2 »THR2 
506  FORMATE5IA3 
GO  TO  507 

505  IF(lGOPT-l)  A20,021  ,522 
*1  22  WRlTt(61 ,523) 

A23  FORMAT ( 1 HI * ASH  THE  PROGRAM  SHOULD  HAVE  C3“I:U  THE  ENTIRE  TAPE) 

GO  TO  A3 r 

A 20  WRITEC61 ,A2A ) 

A2A  FORMAT ( 1H1 . 70HTHE  PROGRAM  SHOULD  HAVE  CO»i:D  ALL  THE  DATE  POINTS 
1WH1CH  FALL  BETWEEN) 

WRITEE61.0253  LATMIN.LATMAX 

A2E  FORMATE  1 Hu , 7X.16H SOUTH  LATITUDE  , FI 0 . 5 , 5 X 16 HNOR TH  LATITUDE  .FIT. 

15  3 

WR1TEE61 ,A26)LONGMIN,LONGMAX 

A 2 6 FORMATE  1 Hj , 7X1 6HWEST  LONGITUDE  , F lU . 5 . 5 X 1 .HE  AS T LON&ITUOC  tclR.5 
13 

GO  TO  A3r 

A 2 1 WRITEC61.A273 

A 2 7 FORMATE 1H1 ,V5HTHE  PROGRAM  SHOULD  HAVE  C0XI2D  ALL  THE  DATA  POINTS  W 
1HICH  WERE  TAKEN  BlTUEEN  THE  SPELIFIEO  0ATt>3 
4 3t  IFEINEEO-235Ci,502,5Al 
501  E NDF ILF  a, 

REWIND  06 
■>f  2 STOP 
END 
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